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ABSTRACT 


The scope of this thesis is to evaluate the use of the MINI RANGER, 
LORAN (TD) and LORAN (DISPLAY) navigation systems in order to support 
the collection of current profiles by an Acoustic Doppler Current Profiler Re- 
corder (ADCP). 

A theoretical error analysis of these systems is undertaken in order to estab- 
lish the minimum error limits which might be expected when averaging current 
profiles over time frames of up to 30 minutes. a | 

Experimental data was collected with all of these systems in the Monterey 
Bay and was analysed, the results being presented in this thesis. In addition GPS 
data was also collected but time has prohibited its analysis and subsequent 
inclusion. 

The results show that because of ship fluctuations in course and speed there 
IS no Statistically significant difference between the navigation systems when we 
average the data over times of about 25 - 30 minutes. However, they also show 
that both the MINI RANGER and Bottom Tracking with the ADCP can 
produce reasonable results in as little as three minutes, although the ADCP re- 


sults are clearly biased. 
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I. INTRODUCTION 


The relative accuracy of navigation data is very important in the processing 

of Acoustic Doppler Current Profiler (ADCP) data. 
_ The shipboard ADCP is an instrument which measures vertical profiles of 
currents. The raw ADCP data is collected in 128 four meter bins, and is used in 
conjuction with a navigation system from which the position of the ship is re- 
corded. 

In this experiment we will try to evaluate which of the navigation systems 
available to us 1s better for use with the ADCP in order to determine the absolute 
Gliemt velocity. We collected data with MINI RANGER, LORAN (TD), 
Bete. (DISPLAY), and GPS receivers. Due to time constraints, however, the 
GPS data has not been included in the analysis. 

Chapter 2 describes the cruise from which the data analysed in this thesis 1s 
taken. It also describes the time svnchonization problems which existed between 
the various tvpes of navigational equipment used. 

Chapter 3 describes the ADCP and its principles of operation. 

im chapters + and 5 we*try to define the theoretical random error for the 
LORAN (LC 408) and the MINI RANGER respectively. 

Chapter 6 gives the results from the data collected in a cruise of Monterey 
Bav while Chapter 7 gives a summary of the conclusions from this study together 


with recommendations for future work. 


I. DATA COLLECTION 


A. CRUISE 

The day chosen for the cruise was 22 September 1989 with the ship R/V Point 
Sur. During the cruise the following equipment belonging to the NPGS and R/V 
Point Sur was used: 

1. Motorola Mini Ranger, 

2. LORAN C receiver (Internav LC 408), 

a= rin pierGie Suncceivicrw 

4. HP Computer logging three data sources simultaneously, and 

5. RDI ADCP with the Data Acquisition System (DAS) and IBM XT com- 
patible computer. 

The results reported in later sections will concern only the MINI RANGER 
and LORAN C navigational systems. 

The ship navigated the runs from point | to point 2, 2-1, 1 - 3,3 - 2, 2 - 3, 
and 3 - 4 (see Figure 1). The speed during the measurements was 6 knots which, 
from previous experience, has proven effective fon ADCP data collecnonm stn 
running from | - 2 for example, a turn was made at point 2 and then the ship 
Steamed back to point 1 without attempting to follow the identical outward track. 

The average depth in the area of the triangle | 2 3 is 30 fathoms, and for the 


run 3 - 4 it varies between 20 to 300 fathoms. 


B. TIME SYNCHRONIZATION 

During the cruise the four pieces of equipment used each had clocks which 
operated independently of each other. In order to synchronize all these clocks, at 
five times during the cruise time comparisons were made between them by using 
the Mini Ranger as the common time base (see for example Table 1). 

These times were then converted from a MR time base to an ADCP time 
base (see Table 2). 

The standard error in the time comparisons was at least 1 second. Once a 


dedicated multi-channel data logger is installed which allows all equipment out- 


Monterey Bay Vicinity 


w2°w 





Figure 1. CHART OF THE WORK AREA 


puts to be recorded on a common time base, it is anticipated that this synchroni- 


zation error Will be reduced to negligible levels. 


Table 1. COMPARISON OF CLOCKS WITH THE MINI RANGER (MR) 


LC 408 15:14:00 15:14:00.6 
ADCP | ye ee ce) 15:15:06.6 










Table 2. DIFFERENT TIMES BETWEEN THE CLOCKS WITH COMMON 
TIME BASE (ADCP) 


LC 408 15:14:00 ADCP 15:14:04.0 


With five tables on a common ADCP time base, a linear regression between 










the different times was performed for each equipment type in order to calculate 


the coefficients A and B of the equation: 
Y=AX+ 8B (1) 


The linear clock drifts were thus calculated and used to assist in interpolating 
the data into the ADCP time base. 


Il. ACOUSTIC DOPPLER CURRENT PROFILERS 


A. INTRODUCTION 
The shipboard Acoustic Doppler Current Profiler (ADCP) from RD Instru- 
ments is a State of the art instrument which measures vertical profiles of currents 
of the ocean at points along the ship’s path [RD Instruments, 1989]. It allows the 
ocean to be sampled in a way which is fundamentally different from moored 
currents or drifters (instruments most commonly used for current measurement). 
Consider the general case of measuring the current at the water parce! whose 


coordinates are X,,(¢) using an instrument whose location is X(t) 
A tila Ap( ar (2) 
Where r is a vector between the two locations and 
ay Cee (f) (3) 
From (2) and (3) 
A(t) = Al) + LAO) — Aol)! (4) 
The current U at the measurement point X, + r will be: 


_ a 
Ai RA eS eek) 5 
ill ella cami 6, salaliaanas F (9) 


Equation (5) can be rewritten as 





4 aXo x 
O(N a) — + V(r) (6) 
dat 
Where V(r) is the velocity of the water parcel relative to the instrument posi- 
tion and dXo.dt is the velocity of the instrument itself with respect to the Earth. 
Since the currents are a difference between two directly measured quantities 


dxo dt and V. the measurement of currents to an accuracy of | cm sec from a 


tn 


ship travelling at 10 knots (approximately 5.2 m/see or 520 em/see) requires that 


both dXo/dt and V must be measured to an aeeuraey of 0.2%. 


B. PRINCIPLES OF OPERATION OF THE ACOUSTIC DOPPLER 
CURRENT PROFILER 

ADCP’s use the Doppler effect by transmitting a succession of aeoustie pulscs 
at a fixed frequeney and listening to the resulting baekseattered water mass ech- 
oes in aS many as 128 depth cells (bins) over a depth range of 30 to 700 meters. 
When seatterers move toward the ADCP, the echoes heard by the seatterers 1s 
Doppler shifted to a higher frequency. The amount of this shift is proportional 
to the relative veloeity between the ADCP (ship) and seatterer Figure 2. 
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Figure 2, WOLUME SCATTERING OF SOUND 


Part of these Doppler-shifted echoes refleet backwards or are backscattered 
to the ADCP. The backscattered echoes appear to the ADCP as if the scatterers 
were the souree of the eehoes (see Figure 3). 

The ADCP hears the baekseattered echoes, Doppler shifted a second time. 
Since the ADCP both transmits and receives, the Doppler shift is doubled, and 
the equation for the Doppler shift will be 

Fy = 2Fs~—e0s A (7) 
where 


e /, 1s the Doppler shift, 
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Figure 3. FIRST-SECOND SOUND DOPPLER SHIFT 


e F, is the frequency of the sound when everything is still, 

e A is the angle between the acoustic beain and the water velocity, 

e V is the relative velocity between the sound source and the sound receiver, 
and 

e Cis the speed of sound (m/sec) in the water at the face of the transducer, 
and can be calculated from the expression 

a ramen |0) i 2910" T " 
+ (1.34 — 1072 7)(5 — 35) + 1.61077D . 
where 


e D = Depth, in meters, 
e S = Salinity, in practical salinity units, and 


e T = Temperature, in degrees Celsius. 


The water temperature at the face of the transducer is measured by the 
ADEPY The speed of sound (C) can be calculated from this measured water 


temperature using mommna values for depth (D) and salinity (S). 
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Figure 4. WATER VELOCITY VECTOR 


Computer analysis of the Doppler frequency shift of backscattered echoes 
from each bin is used to generate a precise depth segmented picture (or profilc) 
of water currents throughout the water column bounded by the path of the 


acoustic beams. 


C. THREE DIMENSIONAL CURRENT VELOCITY VECTORS 

The ADCP beams each measure a single velocity component, 1.e. the compo- 
nent of velocity toward or away from the transduser. When using multiple 
beams, one must make an assumption that currents are the same (homogeneous) 
over layers of constant depth. When the ADCP uses multiple beams pointed in 


different directions, it senses different velocity components. 
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Figure 5. VELOCITY RESOLUTION WITH FOUR ADCP BEAMS 


Figure 5 shows how the ADCP, using four acoustic beams, computes three 
velocity components. The first pair of beams produces one horizontal and one 
vertical component, the second pair produces one horizontal (perpendicular to the 
first horizontal) and again one vertical component. 

The product SCALE x (V, + V, -— V,;— V,) is defined as the error velocity 
component, where I’, to VY, are the velocity components along the beam directions 
Figure 6. This sum should be close to zero. The error velocity component allows 
us to evaluate if the assumption for the horizontal homogenuity is reasonable. 


Also this is another way to estimate data quality. 
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Figure 6. SMALL AND LARGE ERROR VELOCITY 


In Figure 6 there are two different situations, in the first the velocity is the 
same in all four beams at a constant level, in the second the velocity in one beam 
is different. The error velocity in the second case will be on average larger than 
the error velocity in the first case. The difference in the error velocities can be 


from different currents or from errors caused by malfunctioning equipment. 


D. DEPLOYMENT AREA 
The RD Instuments ADCP remotely measures water flow velocity along the 
lines of position defined by the four narrow vertically inclined acoustic beams. 


To insure accurate current measurement it is necessary that the water mass in the 


region of measurement be free of strong acoustic reflecting objects (e.g platform 
members, large cables, surface, etc.) within a +15 degrees conical sector along the 
direction of each of the beams. However, since only three beams are required for 
computation of three axes current components, in applications where potential 
interference objects may be close to one beam, the current vector may be calcu- 
lated from the other three beams. 

In shallow water vertical profiling applications, acoustic interference from the 
Surface (upward looking ADCP) or bottom (downward looking ADCP) limit the 


vertical current measurement region to a maximum range defined by 
Aas — eae Osae (9) 


where 
e R(max) = maximum profiling range, 
e D = distance to surface’ bottom boundary, and 


® = acoustic beam angle relative to the vertical. 


E. BOTTOM TRACKING OPERATION 

The ADCP can be used in a bottom tracking mode to give direct estimates 
of the velocity components of the ship. In the experimental data collected for this 
thesis we used the bottom tracking on for the runs 3 - 2, 2 - 3, and at the first part 
of the run 3 - 4. In order for the ADCP to gather bottom tracking data, bottom 
echoes must be distinguished from other echoes. The ADCP transmits a dedi- 
cated bottom track ping between current profiling pings. The number of profiling 
pings can be Selected by the user. Consequently, the bottom tracking ping data 
is Separated from the current profiling data. Bottom echoes are identified by vir- 
tue of their greater echo strength (the bottom returns a stronger echo than the 
echo from a profiling bin). 

Ship velocity (relative to the bottom) is measured in much the same wav as 
current velocity (the Doppler shift of the backscattered bottom echo being pro- 
portional to the ship velocity). Bottom depth is detected by comparing received 
echo amplitude to a detection threshold which decreases with increasing range to 


the bottom. 


eveeociry PROFILES 

The most important feature of ADCP’s is their ability to measure current 
profiles. ADCP’s break up the velocity profile into uniform segments called depth 
ells. 

Each depth ceil is similar to a single current meter, but there are two basic 
differences. The first difference is that the depth cells in an ADCP are always 
uniformly spaced while current meters can be spaced irregularly, the second 1s 
that the ADCP measures average velocity over the depth range of each depth ccll 


While the current meter measures current only at the current meter Figure 7. 
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Figure 7, ADCP AND CURRENT METERS 


The advantage of the ADCP averaging the velocity over the full range of a 
depth cell is that it reduces the effects of spatial aliasing (Signals at frequencies 


higher than the time series can resolve are mistaken for low frequency signals). 


G. ACOUSTIC DOPPLER CURRENT PROFILER DATA PROCESSING 
The raw data is collected in 128 four meter bins over a three inmute ensemble 
interval. Generally speaking, either a MINI RANGER, LORAN C (calculated 
@ , 4 from time differences or ¢@ , 4 recorded from the display) or a GPS receiver 
are used for navigation and a position recorded every 30 seconds. By interpo- 


lation a position can be calculated at the end of each ensemble. 


The first steps in processing are the calculation of ship’s velocity from the 
navigation data. From this navigation data the U and V components of ship’s 
velocity are calculated. 

The next step in processing 1s the initial determination of the depth to which 
the data of each ensemble remains reliables For our case the, 75 Dl 
chosen. The basic criterion comes from the good percent of return echoes (BIN 
STATISTICS FILE). By subtracting the ship’s velocity from the average velocity 
within the chosen reference layer an absolute reference layer velocity for each 
ensemble is obtained. The series of absolute reference velocities is then filtered 
with a low pass Hamming window filter. 

Once the absolute reference velocity is determined the velocity profiles of each 
ensemble with respect to the reference velocity are also determined, thus vielding 
the final profiles of absolute water velocity. The remaining profiles of absolute 
velocity are then averaged over the time interval for each run. 

In appendix E there are 18 average current profiles from different navigation 


data. 


IV. THEORY OF LORAN C 


A. INTRODUCTION 

LORAN C (LONG RANGE NAVIGATION) is a pulsed, low frequency 
(100 KHZ carrier) long-range hyperbolic navigation system. It operates on the 
principle that the difference in time of arrival of signals from two stations, ob- 
served at a point in the coverage area, is a measure of the difference in distance 
from the point of observation to each of the stations. 

Measurements of the ship’s latitude and longitude ( @ , 4) and Time Differ- 
ences (TD) from a LC 408 LORAN C receiver were recorded throughout the 


experiment (22 September 1989) at intervals of 30 seconds. 


B. SYSTEM DESCRIPTION 

LORAN C_ normally revuires signals from at least a Master and two Sec- 
ondary Stations to give a positional fix. 

Transmitters are grouped to form a “chain” of which one station is labelled 
the Nlaster (M) and the others are called secondaries (X,Y,W) as shown in Figure 
8. 

The chain designators for LORAN C are 4 digit numbers which indicate the 
pulse Group Repetition Interval (GRI) in tens of micro-seconds. For example, the 
west coast LORAN C chain (Figure 8) is designated 9940 and has a GRI 99400 


micro-seconds ( see Figure 9) [LC 408, Operation Manual] 


C. PULSES-PHASE AND CODES-CYCLE SELECTION 
Each station of a LORAN C chain transmits groups of pulses (Table 3). 


Table 3. LORAN C PHASE CODES 
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Figure 8. LOCATION OF WEST COAST LORAN C STATIONS 
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Figure 9. LORAN C GROUP REPETITION INTERVAL (GRIT) AND TIME 
DIFFERENCES (TD) 


The Secondary transmits eight pulses and the Master nine (the last pulse for 
indentification and blink alarm). The pulses are phase coded to improve the sig- 
nal to noise ratio through compression and to distinguish between Master and 
Secondary transmissions. 

Low frequency radio signals propagate over the earth’s surface at nearly the 
velocity of light in a vacuum. However, Maxwell’s equations dictate that the 
ground wave (surface wave) will be influenced slightly by the surface parametci's 
of geometry and electrical properties. 

In order to make the received 100 KHz signal more stable and reliable within 
a given coverage area, the LORAN C radio navigation system is designed as a 
pulse system which separates the ground wave from the skywave. 

Because the earth parameters remain nearly constant, LORAN C has dem- 
Onstrated a repeatability of quite high accuracy. 

The high accuracy of LORAN C, despite long ranges from transmitters, is 
due to a technique called “cycle matching” (Bowditch, 1984]. 

The LC 408 tracks the third cycle cross-over path point Which is very con- 
sistent between transmissions and less susceptible to skywave interference than 
later cycles. By tracking this cycle on all pulses, high accuracy is attained. Figure 


10 shows the third cycle tracking point with an example of skywave interference. 
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Figure 10. THIRD CYCLE TRACKING POINT 


D. TIME DIFFERENCE EQUATIONS 

The elapsed time between the arrival of the signals from the Master and a 
Seeondary stations 1s called the “Time Differenee” (TD). Eaeh observed time 
differenee, or rate, provides one hyperbolie line of position. By observing the 
transmissions from four stations, three hyperbolie rates are measured and the 
position ean be determined by either graphical or analytical techniques. 

The time difference observed at a receiver is the difference in arrival times for 
Signals from the master and one secondary transmitter in the ehain. Beeause all 
transmitters share the same frequencies, their signals must be separated in time 
to prevent interference. 

The chain is syehronized so that the Master transmits first followed by each 
of the Secondaries. The transmission of each secondary is specified by the emis- 
sion delay so that in the coverage area signals from one station will overlap an- 
other [Schenebele, 1979]. 

Suppose that at a point P there is a reeeiver. The observed time difference 


using the Master M and Seeondaries X,Y,W are 
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Figure 11. TYPICAL HYPERBOLIC FIX GEOMETRY 
HO; 1h” = EDyy Ae Cy , {a4 
yg by = ED, + ly aoe lyf 
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where ED,, EDy, EDy are the emission delay for W, Ay Y. The termsai sees 


are travel times from the W , X , Y to P. The term ¢, is the travel time from the 
Master M to P. 
In order to express the time difference as a function of geographic position , 


the travel time t is Separated into additive terms 


Xx D+F (13) 


WNehe: 
e C = free space propagation velocity, 
¢ n = index for refraction for a standard atmosphere, 
e D = geodetic distance from the transmitter to receiver, and 
e F = phase factor which corrects for effects of the earth’s surface along the 


path. 


Substituting the last formula into the equations (10), (11), (12). gives the 


equations: 
TDi =o a (Dye = ID) = Fe (14) 
[DX = ED a x (Dy — Dy) + Fy — Fy (15) 
TDY = ED eee (Dy eee ane (16) 


All these equations relate the time differences WD (m4 y to hema 
tances from each of the four transmitters. 

Since the latitude and longitude is mathematically calculated from the time 
difference numbers, the crossing angles, gradients, and signal to noise ratios are 


Still important. 


E. HYPERBOLIC GRADIENT, CROSSING ANGLES 
For a given pair Master, Slave a family of hyperbolas separated by a costant 


value can be plotted. 
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Figure 12, HYPERBOLIC GRADIENTS 


The hyperbolas in Figures 12 and 13 are all separated by 10 microseconds 
[Gann]. That means that the distance on the baseline represents about 1852 
meters. The distance between successive hyperbolas increases as one moves io- 
wards the baseline extensions or away from the baseline. This change in the ac- 
curacy of a hyperbolic Line of Position (LOP) that occurs relative to the position 
of the LORAN C stations is known as the gradient. 

The angle between the LOPs 1s called the crossing angle Figure 13. The 
crossing angle must be bigger than 30 degrees and smaller than 150 degrees. T luis 
fact is common with any ranging positioning system, but is further limited for the 
LORAN C because the reception of three stations is required to define two LOPs 


(direct ranging navigation systems provide one LOP per station). 
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Figure 13. HYPERBOLIC GRADIENTS AND CROSSING ANGLES 


F. ERROR ANALYSIS FOR LORAN C 

Any discussion of accuracy should begin with the definition of the types of 
accuracy. 

Repeatable accuracy (Repeatability, Precision) is the accuracy with which a 
navigator can return to a location, the coordinates of which have been previously 
measured using the same system. 

Random error is one which results from basic limitations in the method. The 
characteristics of this error can be determined by statistical analysis of a sufficient 
number of measurements. This type of error affects the repeatability. These cr- 
rors are identified by their Gaussian distribution. 

Peele accuracy is the accuracy with which a navigator can determine his 


point position in terms of a Well defined coordinate system e.g.(@ , 2 ). 
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Systematic errors result from a basic (but unrealized) fault in the method 
and cause the values to be consistently biased from the true value. It affects the 
absolute accuracy, and often cannot be detected by statistical analysis. 

The LC 408 microprocessor uses an exacting spherical trigonometric formula 
to derive its latitude and longitude positions. Also it is important to note that the 
LC 408 used for the calculations assumes that the LORAWN signals travel over 
ideal all seawater paths from the transmitter to the receiver. Tp correct for path 
over land, the LC 408 allows use of “additional secondary factors”; this is typi- 
cally used only in small regions adjacent to known fixed locations. 

For the error analysis given shortly, the shape of the earth is approximated 
by a sphere, and it 1s assumed that all the west longitudes are aan. 

In order to find the standard deviations in the velocities, we established a 
fixed point in the roof of Spanagel Hall to enable an assessement of the standard 
deviations in the three time differences and after it to propagate the error at the 
positions and the velocities. 


Consider the spherical triangle PSM (see Figure 14 [Cross, 1981]). 


> 


Figure 14. SPHERICAL TRIANGLE SPANAGEL-M-NORTH POLE 
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Where (M) can be the Master, S can be the point on the Spanagel roof, and 
P can be the North pole. We try to find an expression for the distance between 
Master and Spanagel roof with respect to @ , 4 (coordinates of the point on the 


Spanagel roof), 


m=90—-—¢ (17) 
s=90- Py (18) 
P=)—)y (19) 


Where @,, and Ay are the coordinates of the Master station. Using the law of 


cosines for the side p, 
cos p = cos mcos s + sin msin scos P (20) 
Substitute (7) , (18) 719) ites 20) io sceia 
cos p= sin Py, sin @ + cos Py, Cos @ cos(Z — 44,7) (21) 
Differentiate each part of (21) 


— sin pdp = cos $y;sin ddhy, + sin hy, cos Pde 
— Sin @y, cos @ cos(s — Ay )dPyy 
= COS @,,Sin © cos(7A — Age (22) 
— cos dy,cos msin(d — Ays)dZ 


“aaeeOs Oy COs “Dp sin(/ = A)d/ x4 


Since the position of Master station is a fixed point d/, = doy = 0 equation (22 


becomes 


— sin pdp =[ sin @4, Cos @ — cos @y;sin Pcos(4 — Ay,)|d@ 


= | (23) 
— cos Py cos Psin(A — Ay, )d/ 


We now introduce into each term in the right part of (23) the factor sin p. Using 


the relationship between two angles and three sides 
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SiycOSS — "COs 5 Sion — Sin seas? COS (24) 


Substituting for s, m and P as appropriate, we have 
Sin Py cos @ — cos day sin d cos(A — Ay) = 
cos(90 — #44) sin(90 — ) — sin(90 — dy) cos(90 — @) cos(A — Ay) = (25) 
sin pcos S$ 


From the law of sines we know that 


sinm — Sin p 


sin.if sin P oO) 
OR 
Sin 77sin P =sin psin M/ (28) 
Substituting for m and P and multiplying both sides by cos dy , we have 
COS Py, COS @ Sin(4 — 447) = 
COS G1, Sin(90 — o)sin(/ — 4,,) = (29) 
sin p sin Af cos }yy 
Substituting (25) and (29) into (23), we have 
dp = — cos Sd + sin Af cos @yd/ (30) 


Actually in order to convert the angle dp into distance da (arc length), using 


Spherical approximation multiply by the mean radius of the earth R 
da = — Rceos Sd@ + Ksin M cos Pyjd/ (31) 


From the last equation we have 


: | 
=, = — Roos $ (32) 


ty 
ow 


Ca 


nee Ros @y sin M 


(33) 


where M and §S are the interior angles of the triangle PSM. In similar way can 


form the equations for the other spherical triangles 
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The sign changes in the second equation due to different geometry 
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Next let 
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Suppose that the velocity of the the electromagnetic propagation is constant and 


eq @) them: 


b> 
I 
|S 


—_— 


ad, 
ac 
ds 
on van 


(43) 


(44) 


(45) 


We propagate the errors into these expressions by forming the Jacobian ma- 


trix. We assume here that time is subject to negligible error. This matrix is given 


bv 


and 
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(47) 


(49) 


od, Cal, A) 6a,(@, A) 


a (a 
= Reos dysin M+ Rceos @ysin X (50) 
= R( cos @y sin M + cos dy sin X) 


Gds  ald.4) ~~ Ba(, A) 


mec | es 
= Rceos dysin M — Reos dysin Y (51) 


= R( cos dy sin M —cos ¢gysin Y) 


bd; dal$, 4) Gax(, 4) 


nn) Sny, 
= Rceos Py sin Af — Reos my sin IV (52) 
= R( cos @y sin Mf — cos Py sin IV) 


C=3x 10°(A/ sec) = 3 x 107 '(.M/nsec) 


Now define the variance covariance matrix of LORAN time differences, as- 


suming that there is no correlation between the three time differences. 
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The LC 408 is provided with control of the filtering time constant. Selectable 
time constants are 0 seconds for no filtering, 5 seconds for good response, 10 
seconds for slower response, and 20 seconds for slow craft or monitor applica- 
tions. 

In order to measure the standard deviations for the time differences, we re- 
corded data for four days on Spanagel Roof, each day (24 hours) with a different 
time constant. This resulted in four different variance covariance matrices. 
There are all given in Appendix A. 

From these results we conclude that the filters dont work for the time vai 
ference between Master and Secondary Whisky because the standard deviation 
remains approximately the same for the different time filters. We have no ex- 


planation for this beyond it being an instrumental problem. We also conclude 


that for monitoring applications, the longer the filtering time, the smaller the 
Standard deviations. 


ive weight matrixeP = >-! is given by 


l/c? 0 0 
P=| 0 I/o° 0 (55) 
0 0 1/o?, 


The variance covariance matrix for the position [Uotila, 1986] will be: 
Ep, y= (Ux Px J (56) 


iWemme the matrix Dp) as 


S=15 (37) 


Now we form the biock matrix &£& ( $,/,@,/,) which is the variance 
coVarlance matrix for two positions (the error characteristcs are the same for the 
two positions). If 1t 1s assumed that there 1s no correlation between the two posi- 


tions, then 


ee Or 5 
S64 /ids/2)=| = v (4. (58) 


Ay) 


In order to convert the difference in latitude and difference in longitude (radians) 


between two positions in meters, apply 
ANG XK (59) 
BaeN/ Xo COSs ree ck (60) 


where @,, is the mean latitude of the work area, and R 1s the mean radius of the 
earn in) 6371 KM. Let 


C,=cos ¢,, x R (62) 


The components of the ship velocity V , U (Northing, Easting) are then given by: 


c 
Ve == x (2 $1) (63) 
C 


The Jacobian matrix for the velocity will be 
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Consider time intervals between positions 3 ,6,12, 20 and 40 minutes. The 


Variance covariance matrix for the velocity components will be: 


( 
v ae bi 
op y=d x ey ree es yi (67) 


Appendix B gives 5 tables for standard deviations and covariances of velocity 
components with different time intervals between positions. From these tables 
we see clearly that as the time interval between the positions increases the stand- 
ard deviations for the V and U components decreases (assuming constant course 
and speed of the ship for the given interval). Also we note that when time aver- 
aging over 20 minutes or longer, the choice of filter is not critical. 

Appendix C gives a program (DRIVLR FORTRAN) written by the NOAA 
Which calculates the geographic position of the ship from two LORAN Time 
Differences (TD) and a Dead Reckoning (DR) position of the ship (inverse com- 


putation). It can also be used to convert the geographic position of the ship into 
LORAN Time Differences (forward computation). 


V. THEORY OF MINI RANGER 


A. INTRODUCTION 

The Mini Ranger Falcon was used during the cruise for navigational posi- 
tioning. The standard Mini Ranger operates at up to 37 kilometers (about 20 
nautical miles) with a probable range measurement error of 2 meters. 

It operates at microwave frequencies and requires that line of sight be main- 
tained between the reference stations and the receiver transmitter. Significant 
obstuctions such as land masses, buildings, or dense foliage will interfere the op- 


eration of the system. 


B. REFERENCE STATIONS 
The positions of the reference station sites in the UTM coordinate system are 
[Krioneritis. 1989] listed in Table 4. 


Table 4. REFERENCE STATION SITES 
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C. RANGE POSITIONING 


The Mini Ranger determines a two range position when the lengths of all 
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three sides of a triangle are known (trilateration). Suppose that a certain task re- 
quires activity between points © , DE (Pigiire 1S )sthe relencmee sonics nen am 
points A and B. 1n the acceptable area the angle of intersection between the two 


range lines is between 30 and 150 degrees. 
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Figure 15. PLANNING THE MOST DESIRABLE SITE 


Arcs AXB and AYB define the maximum and minimum distance to maintain 
between 30 and 150 degrees respectively. 

Generally the Mini Ranger uses the method of least squares to calculate po- 
sitions from three or more reference stations. Operation in areas where the ge- 


ometry is poor will degrade positional accuracy [Mini Ranger operation manual, 


1981]. 


D. ERROR ANALYSIS FOR MINI RANGER 
While the Mini Ranger measures ranges well, the actual position error re- 


sulting from a set of range readings depends upon a number of contributing fac- 


Ole: 


1. syStem error, and 


2. geometry changes. 
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Figure 16 shows the maximum positional error for two range geometries of 30 , 


90 and 150 degrees [Mini Ranger operation manual, 1981]. 


STANDARD ERROR 


POSITION ERROR = 
sin @ 


(68) 


The position error for various two range geometries, assuming range error of a 


single measurement of 2 meters are given in Table 5. 


Table 5. POSITION ERROR FOR DIFFERENT 
CROSS ANGLES 












For this experiment we'll take the information for the variance covariance 
matrix of the positions from the Mini Ranger data processing (see, for example, 


Krioneritis.1989).. Deliment nematic aas 


0 0 
S=|9 of (69) 
The variance covariance matrix for the position will be 
>( 6 y) 6? Ce “ie (70) 
wet HN gy nag OV y Gy 


We now form the block matrix 2 ( 4, ¥,X;,Y;,) which is the variance 
covariance matrix for two positions. If it is assumed there is no correlation be- 


tween the two positions, then 


at 2 meters 


REFERENCE 


nrtatn ten 





Figure 16. MAXIMUM POSITION ERROR 


ie £(A), ¥) 20 : 
LX, } Sh: >) — a X(X5, 5) ( { 1) 


Define the velocity components as 
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As we did with LORAN C, we propagate errors into these expressions by forming 
the jacobian matrix for the velocity components. We again assume that time is 


subject to negligible error. This matrix is given by 
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=<; [1 ee 1 bo 7) 
The variance covariance matrix for the velocity components will be 
if 

ay U = of x 2y, ¥, (Xd, ¥5 x y (77) 


Table 6. MIAX-MIN VALUES OF STANDARD DEVIATIONS AND 
COVARIANCES FOR VELOCITY COMPONENTS FOR 3 MINUTE 
TIME INTERVALS BETWEEN POSTTIONS (MINI RANGER) 
MAX-MIN 


Table (6) shows that the the standard error in the velocity components 







doesn't vary greatly given the geometry 1n the work area. The results come from 
6 different points in the whole work area. 
Appendix D gives the velocity error results for different time intervals. Simi- 


lar to the LORAN (Appendix B), the standard deviations for the velocity com- 


ponents decrease when the time interval between the positions increases 
(assuming constant course and speed of the ship for the given time interval). 
Comparing the theoretical results from Appendices B and D, we can see that 
the MINI RANGER is not only a more stable navigation system than the 
LORAN C, but also that it provides in 3 minutes the level of accuracy in the U 
and V velocity components which require 20 minutes to obtain with LORAN C. 
Because of the practical difficulties in maintaining a constant ship’s heading and 
speed for 20 minutes, and thus the need for an averaging process over Such a time 
interval, the value of using a MINI RANGER wherever possible is even further 


enhanced. 


oh 
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VI. RESULTS 


In this chapter there are eleven Tables which summarise the results of the 
computation of the ship velocity components from MINI RANGER, LORAN 
Time Differences (TD), LORAN (DISPLAY) and Bottom Tracking. We choose 
to separate the LORAN (TD), and LORAN (DISPLAY) results because of the 
rounding off which occurs when using the LORAN (DISPLAY). 

The LORAN results given here were obtained using different filters (2 - 1, | 
- 3 filter with time constant 5 seconds, | - 2 filter with time constant 20 seconds, 
3 - 2,2 - 3 filter with time constant 10 seconds, 3 - 4 filter with time constant 0 


seconds ). 


Table 7. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAY) FROM POINT 1 TO POINT 2 AT 3 MIN- 
UTE TIME INTERVAL. 













LORAN (DIS- 
PLAY) 


TIME 


MINI RANGER EORAN GG) 


36 


Table 8. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAY) FROM POINT 2 TO POINT 1 AT 3 MIN- 
UTE TIME INTERVAL 













LORAN (DIS- 


MINI RANGER LORAN 
oe PLAY) 
TIME 
14 57 10 | 265.3 | 134.3 | 254.5 | | 135.0 | 258.0 
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feo | vse | ae | wr [aes [in | a 
CER ss [we [oes [sr [aa 


From Tables 7 and 8 we can sce a discontinuitv beginning at 14 36 11 near 





Py 





the end of leg 1 - 2 and ending at 15 03 10 at the beginning of Ieg 2-1. This 
discontinuity occurs at the same position down each leg and 1s of approximately 
the same magnitude in each case. While it is tempting to view this discontinuity 
as being due to a surface current, the inconsistency in sign prohibits such an ex- 
planation. It can only be concluded therefore that it is due to a ship navigation 
correction. 

In these tables we see clearly that the LORAN (DISPLAY) results show a 
much geater variability than the LORAN (TD). Also if we calculate the mean 
value and the standard deviation of the mean of the absolute difference of the U 
and V components between the MINI RANGER and LORAN (TD), and the 
MINI RANGER and LORAN (DISPLAY). we have Table 9. 


Table 9. COMPARISON BETWEEN MINI RANGER AND LORAN (TD) AND 
THE MINI RANGER AND LORAN (DISPLAY) FOR LEGS 1 - 2 AND 


2- i. 









MIINT RANGER - LORAN (TD) MINI RANGER -LORAN (DISPLAY) 


LEG 
AU (¢,y) AV (G,y) A U (¢yy) AV (6,y) 
2.9(0.9) 4.5(2.1) 2.9(0.9) 7.9(1.8) 
1351053) 4.5(1.2) 4.2(1.6) oh) 






After using the FISHER-BEHERENS test [Hamilton, 1964] to compare AU from 
MINI RANGER-LORAN (TD) against AU from MINI RANGER-LORAN 
(DISPLAY) and AV against AV we see that these are statistically different at 
95% confidence interval in 3 out of 4 cases. This provides strong evidence to 
Suggest that the results from LORAN (TD) and LORAN (DISPLAY) are from 


different statstical populations. 


Table 10. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAY) FROM POINT 1 TO POINT 3 AT 3 
MINUTE TIME INTERVAL 









LORAN (DIS- 
PLAY) 





MINI RANGER LORA NTIS) 















U V g V U V 
(cm.sec) | (cm.Sec) | (cm/sec) | (emisec) | (emmsec)) |i emmeae 
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Table 11. COMPARISON BETWEEN MINI RANGER AND LORAN (TD) AND 
THE MINI RANGER AND LORAN (DISPLAY) FOR LEG 1 - 3 


NMIIN]I RANGER - LORAN (TD) MINI RANGER -LORAN (DISPLAY) 
LEG 
1.9(0.6) 2.9(0.6) 7.0(1.7) 7.9(1.7) 












Table 10 shows similar results for leg 1 - 3. From Table 11 it is also obvious 
that the results which comes from the LORAN (TD) and the LORAN (DIS- 
fer rare different for leg | - 3. 


Table 12. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAY), BOTTOM TRACKING FROM POINT 3 TO 
POINT 2 AT 3 MINUTE TIME INTERVAL 













LORAN (DIS- BOTTON 


N AN 
MIINT RANGER LORAN (TD) PLAY) TRACKING 
TIME 
161618 | 103.5 | | 103 | gal | 1084 116.5 | 5 
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Table 13. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAY), BOTTOM TRACKING FROM POINT 2 TO 
POINT 3 AT 3 MINUTE TIME INTERVAL 











LORAN (DIS- BOTTONI 
PLAY) TRACKING 


TIME 


fr ec 
ms A 


MINI RANGER LORAN (TD) 


From Tables 12 and [3 we can see that the more consistent results comes 
from MINI RANGER. The variability of the LORAN (TD) results are smaller 
than the variabilitv of the LORAN (DISPLAY) results. We also see that the 
Bottom Tracking, while giving values to the same internal consistency as the 


MINI RANGER, gives results which are systematically biased. 


Table 14. COMPARISON BETWEEN MINI RANGER AND LORAN (TD) AND 
THE MINI RANGER AND LORAN (DISPLAY) FOR LEGS 3 - 2 AND 


MIINT RANGER - LORAN NMIINT RANGER -LORAN MINT RANGER - BOT- 
(TD) (DISPLAY) TONT TRACKING 


LEG 
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From JTable 14, after using the FISHER-BEHERENS test to compare AU 
from MINI RANGER-LORAN (TD) against AU from MINI 
RANGER-LORAN (DISPLAY) and AV against AV we see that these are sta- 


tstically different at a 95% confidence interval in 2 out of 4 cases. 


Table 15. SHIP VELOCITY COMPONENTS FOR MINI RANGER, LORAN 
(TD), LORAN (DISPLAY). BOTTOM TRACKING FROM POINT 3 TO 
POINT 4 AT 3 MINUTE TIME INTERVAL 
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Table 16. COMPARISON BETWEEN MINI RANGER AND LORAN (TD) AND 
TOM TRACKING 


THE MINI RANGER AND LORAN (DISPLAY) FOR LEG 3 - 4. 
MINI RANGER -LORAN 
EEG (TD) (DISPLAY) 
co NG) AN (e,)) AUG.) AV (oa, 
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MINT RANGER - BOT- 





Table 15 again vields similar results for leg 3- 4. From Table 16 it is again 
obvious that the results which comes from the LORAN (TD) and the LORAN 
Meroe y ) arevdifferent. 

In Table 17 are the average velocity components for the different positioning 
systems and for the bottom tracking for each leg. From this table we can see: 


e Data averaged over times of about 18 minutes (minimum) show no statis- 
tically significant difference between MIN] RANGER, LORAN (TD) and 
LORAN (DISPLAY). Any one of these positioning systems can be used for 
these averaged times. This is due to the fact that the variations of the ship 
course and speed in the averaged period of time completely swamp the dif- 
ferent precision which the different navigation systems have in their calcu- 
lation of the U and V components of the ship. The minimum time interval 
from which the averaged data over the different navigation systems become 
the same Statistically depends upon the ship’s variations in course and speed. 
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¢ The results from Bottom Tracking are clearly systematically different from 
those of the others systems. They do, however, show the same level of 
internal consistency as both the MINI RANGER and the LORAN (TD) 
results. 


¢ The much larger standard deviations on the LORAN (DISPLAY) results 
are indicative of the much greater variability in computed U and V compo- 
nents when using this system. 


Table 17. AVERAGE VELOCITY COMPONENTS FOR MINI RANGER, 
LORAN (TD), LORAN (DISPLAY), BOTTOM TRACKING 


LORAN (DIS- BOTTOM TRACK- 
PLAY) ING 


-117.7(3.4)| -279.9(1.2)| -117.2(3.3)} -276.7(2.1)| -116.8(3.1) asooa] =f 
ca 

146.1(3.7) | 256.3(1.9) | 144.7(3.8) | 252.7(1.4) | 143.3(3.2) non] = fo 
-270.1(1.1)} -122.3(1.1)] -268.9(1.4)| -121.2(1.3)] -268.5(2.9) ee ee 


MUNI RANGER LORAN (TD) 


In Appendix E there are 12 figures which show the average current profiles 





from point | - 3,3 - 2, 2 - 3,3 - 4 with the three different navigation data from 
MINI RANGER, LORAWN (FD), LORAN BOISE Paw a itiese plotseswons 
the same thing as the Table !7 1e. no significant difference between the naviga- 
tion svstems when we average over times of about 20 - 30 minutes. 

In Appendix F there are 6 plots for points | - 2, 2 - | but not averaged for 


the Whole length (just plots for specific times between the points). These plots 
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show, as expected, that the U and V components of the current are different for 


different navigation systems. 


VH. CONCLUSIONS AND RECOMMENDATIONS 


From this thesis we can conclude that: 


Because of the variability in the LORAN (DISPLAY) results and the need 
for a reasonable measure of statistical strength, it is not recommended that 
LORAN (DISPLAY) be used for giving averaged results of under 20 min- 
utes. Indeed, it 1s preferable to use LORAN (TD) rather than LORAN 
(DISPLAY) whenever possible. Furthemore, LORAN (DISPLAY) results 
not only show a much greater variability than the LORAN (TD) results. but 
also at a 95% confidence interval consistently fail the statistical hypothesis 
that they come from the same population. 


MINI RANGER results as expected are very consistent. It is felt that good 
ADCP current estimates could be derived using MINI RANGER observa- 
tions for periods as short as 3 minutes. 


When averaging for 20 - 30 minutes all the navigation systems give the same 
RESUS: 


Without averaging we recommended that the MINI RANGER be used 
whenever possible. If this system 1s not available, LORAN (TD) should be 
The Tete molec: 


The theoritical investigation revealed that when time averaging over 20 
minutes or longer. the choice of the filter to be used for the LORAN system 
was not critical. The cruise results appear to support this contention. 


Bottom Tracking results are svstematically biased, but have the same 
Standard deviation as the MINI RANGER results. 
In looking to the future work on this subject it is recommended that: 


Further work is required to understand and resolve the Bottom Tracking 
bias problem. 


If this cruise is done again we suggest checking the LORAN filters by mak- 
ing one long run, dividing it into 30 minute segments, and using a different 
filter on each segment. 
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APPENDIX A. STATISTICAL ANALYSIS FOR THE LC 408 DATA 


The provided numbers are in nanoseconds . 


(1) TIME CONSTANT 0 SECONDS (NO FILTERING) 


STATION N MEAN STANDARD 

DEVIATION 
WHISKY 1440 16312005.09 46. 85663572 
X-RAY 1440 27511618.59 22.27564796 
YANKEE 1440 42735925.76 24.94991412 


MINIMUM MAXIMUM STD ERROR 
VALUE VALUE OF MEAN 


toot / 60,00)  163517160.00. daze 3077 
27511540500 2751172020060 58701467 
42735810.00 42736030.00 0.65748797 


(2) TIME CONSTANT 5 SECONDS (GOOD RESPONSE) 


STATION N MEAN STANDARD 
DEVIATION 
WHISKY 1440 16312033.94 51. 75565604 
X-RAY eee027511616.71 19.-82896815 
YANKEE We60G2735928.43 23. 39276554 


MINIMUM MAXIMUM STD ERROR 
VALUE VALUE OF MEAN 


16311860.00 16312180.00 1.36388129 
ZT Ise OO Me2 / S1V680MGO: 0252253919 
42735860.00 42736020.00 0.61645403 


Peiwere CONSIAN F110 SECONDS {SLOMEK RESPONSE) 


STATION N MEAN STANDARD 
DEVIATION 
WHISKY 1440 16311981.04 49.40674807 
X-RAY Pees 2 7511616. 56 13), (6653231 
YANKEE 1440 42735918.06 18.02041879 
ore ‘CONSTAN E20 
PeererCAlTION S) 
STATION N MEAN STANDARD 
DEVIATION 
WHISKY ial 16311990. 54 46587615596 
X-RAY eM 2751618, 60 12799066592 
YANKEE 1 42735915. 31 1697 4596778 
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MINIMUM MAXIMUM STD ERROR 
VALUE VALUE OF MEAN 


162 iwes0, 00 1631211000 - 30196213 
27511 580R00 ney See S0n00: 066277998 
42735850.00 42735990.00 0.47487973 


SECONDS (\ILONITOR 


MINIMUM MAXIMUM STD ERROR 
VALUE VALUE OF MEAN 


16371800800 16312150700 1.:24528015 
Polis O0 627511660700. 0234510122 
27S 5650, 00 42755980.00 0.44486202 


APPENDIX B. TABLES FOR STANDARD DEVIATIONS AND 
COVARIANCES OF VELOCITY COMPONENTS FOR LORAN C (LC 408) 


Table 18. TIME INTERVAL BETWEEN POSITIONS 3 MINUTES 


Je) INE 
6.8 cmysec -61.184 cm?/ sec? 
Gulmem Sec -51.047 cm? sec? 


10 sec ; aS Ae tCitl, SCC -27.447 cm?/ sec? 
20 sec 4.1 cmysec -24.029 cim?*/ sec? 


Table 19. TIME INTERVAL BETWEEN POSITIONS 6 MINUTES 


Eli 
Se cela Mae 3.4 cmysec -15.296 cin?] Sec? 





5.4 cm,;sSec 3.0 cimnsec =) 2.761 Chi} See: 
4.1 cm. sec 2alnC ii See -6.861 cin?*/ Sec? 
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20 sec Brsecmi/sec 2.0 cm/sec -6.007 cm?/ sec? 





Table 20. TIME INTERVAL BETWEEN POSITIONS 12 MINUTES 


2) cmisec 


Ze lomil Sec 
Ie Cmisece Gre Sec —1.502cn1?/ sec? 


Table 21. TIME INTERVAL BETWEEN POSITIONS 20 MINUTES 


















1.7 cm/sec —3.824cm-*/ sec? 












leoee mi SEC —3.190cr17/ sec? 












Lol eiisec —1.715cm?] sec? 







—1.376cr1*/ sec- 


—1.148cr-/ sec* 
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le2 cm/sec —0.617cm?/ sec? 
i. emysee —0.540cm?/ sec? 





Table 22. TIME INTERVAL BETWEEN POSITIONS 40 MINUTES 


—0.344¢nm?/ sec? 


0.4 cm/sec — 0.287cm?/ sec? 


0.3 cmysec — 0.154¢1717/ sec? 


O73 cmysee — 0.135ci?/ sec? 
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APPENDIX C. PROGRAM DRIVLR FORTRAN 
A. GENERAL REMARKS 


This program calculate the hyperbolic coordinates of a ship expressed in 
LORAN time differences from geographic potisions forward computation , and 
the inverse computation. 

Forward computation: 

Seitput parameters = ALAT, ALON, L, DEL, ISYS, IB, IVEL. 
Sout parameters - 1D, BETA, BLEN, BLEM, IER. 


Inverse computation: 


Sein parameters - ALAT, ALON, L, DEL, TD, ISYSPISENT, ISW, IB, 
EL. 


¢ Output parameters - ALAT, ALON, BETA, BLEN, BLEM, ISENT, IER, 
a 


B. DESCRIPTION OF PARAMETERS: 
ALAT - Latitude array in seconds of arc: 
e ALAT(1I) - Latitude of master station. 
¢ ALAT(2) - Latitude of slave 1 station. 
e ALAT(3) - Latitude of master station. 
e ALAT(4) - Latitude of slave 2 station. 
e ALAT(5) - Latitude of ship position. 
e ALON(1) - Longitude of master station. 
e ALON(2) - Longitude of slave 1 station. 
e ALON(3) - Longitude of master station. 
e ALON(4) - Longitude of slave 2 station. 


e ALON(S) - Longitude of ship position. 


L - Spheroid code array. 
L(1) - Spheroid code array for first rate: 
e Ld) = 1 WGS 1984. 


1) 


° L(l) =2 
¢ L(1) =3 
° Lil) =4 
e Lil) =5 
¢ L(1) = 6 
e Lil) =7 
e Lil) =8 
e Lil) =9 


Bessel spheroid. 
Clarke 1858 spheroid. 
Clarke 1866 spheroid. 
Clarke 1880 spheroid. 
Everest spheroid. 
Fischer spheroid. 
International spheroid. 
WGS 1972 spheroid. 


L(2) - Spheroid code for second rate: 


¢ (2) = 1 
¢ L(2) =2 
e L(2) =3 
e L2)=4 
° L(2) = 
© (2) 
© L(2) 
¢ L(Q2) = 
e L(2) = 


5 
6 
7 
8 
9 


WGS 1984. 

Bessel spheroid. 
Clarke 1858 spheroid. 
Clarke 1866 spheroid. 
Clarke 1880 spheroid. 
Ev crests seinererd: 
Fischer spheroid. 
International spheroid. 
W Gisal?7 2 Sspitcroid: 


DEL - Coding delay array, in microseconds: 


© DEL(1) 
© DEL(2) 


Coding delay for first rate. 


Coding delay for second rate. 


TD - TIme difference array, in microseconds: 


© TD(1) 
© TD(2) 


Time difference for first rate. 


Time difference for second rate. 


ISYS - Loran system array: 


e ISYS(1) 


© ISYS() 


Loran system for first rate. 
ISYS(1) =A: Loran A 
ISYS(1) =c: Loran C 


Loran system for second rate. 
ISYS(2)=A: Loran A 
ISYS(2) =c: Loran C 
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BETA - Baseline delay array: 
¢ BETA(I) — Baseline delay for first rate. 


e BETA(2) _ Baseline array for second rate. 


BLEN - Minimum time difference array, in microseconds: 
e BLEN(1) Minimum time difference for first rate. 


¢ BLEN(2) Minimum time difference for second rate. 


BLEM - Maximum time difference array, in microseconds: 
* BLEM(1) Maximum time difference for first rate. 


e BLEM(2) Maximum time difference for second rate. 


ISENT - Sentinel for DR position: 
e ISENT=0 ~~ DR position not given. 
e ISENT=1 ~~ DR position given. 
ISW Operation code: 
e ISW=!]1 Forward computation. 
°e ISW=2 Inverse computation. 
IB - Switch for baseline computation: 
e IB=0 Baseline computation off. 


e IB=] Baseline computation on. 


ie Ree Error code: 


° IER=1 Invalid velocity of propagation code. 
e ITER=2 Invalid spheroid code. 

e JIER=3 Invalid configuration. 

e IER=4 Invalid operation code. 


e IER=5 Invalid readings. 
e IER=6 Invalid DR position. 


IT - Iteration code. 
IVEL - Velocity of propagation code array: 


IVEL(1) - Velocity of propagation code for first rate: 
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©) €)-O27 42°07 0 


299.6929 kilometers per microsecond 


299.6911624 kilometers per microsecond 


299.708 kilometers per microsecond 


IVEL(2) - VElocity of propagation code for second rate: 


e 0 = 299.6929 kilometers per microsecond 


10 


20 


l 
2 


st 


= 299.6911624 kilometers per microsecond 


= 299.708 kilometers per microsecond 
Tita icc irichriae tricia i noe ree TO arene A Se 
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DESCRIPTION OF PARAMETERS 
DATA DECLARATION 
IMPLICIT REAL*8(A-H,O-Z) 
REAL*8 ALAT(5),ALON(5),DEL(2),TD(2),BETA(2),BLEN(2),BLEM(2),LN 
CHARACTER*1 ISYS(2) 
INTEGER H,M,L(2),IVEL(2),ISYS(2) 


ALAT( 1)=3600. OD0*39. ODO + 33. OD0*60.0D0 + 6. 740D0 
ALAT(2)=3600. ODO*38. ODO + 46. ODO*60. ODO + 57. 110D0 
ALAT(3)=3600. OD0*39. ODO + 33. 0D0*60.0D0 + 6. 740D0 
ALAT(4)=3600. OD0*35.0D0 + 19. O0D0*60.0DO + 18. 3050D0 
ALAT(5)=3600. ODO*36. ODO + 45.0D0*60.0D0 + 0. ODO 
ALON( 1)=3600. ODO*“118. ODO + 49. 0D0*60. ODO + 55. 8160D0 
ALON(2)=3600. OD0*122.0D0 + 29. O0D0*60. ODO + 43.9750D0 
ALON( 3)=3600. OD0*118. ODO + 49. 0D0*60. ODO + 55. 8160D0 
ALON( 4)=3600. ODO0*114. ODO + 48. 0D0*60.0D0O + 16. 8810D0 
ALON(5 )=3600. ODO*121. ODO + 55. 0D0*60.0D0 + 0. ODO 
LGly = a 

L(2) = 1 

DEL(1) = 27000. ODO 

DEL(2) = 40000. ODO 

ISYS(1)=1 

ISYS(2)=1 

ISENT=1 

ISW=2 

IVEL(1)=1 

IVEL(2)=1 

IB=1 


CALL LORAN (ALAT,ALON,L,DEL,TD,ISYS,BETA,BLEN, BLEM, ISENT, ISW,IB, 


‘IER, IT,1VEL) 


DO 200 1=1,914 

READ(7,10) H,M,S,TD(1),TD(2) 

FORMAT( 12,1%,12,1X,65. 2,21x 6.2, 1X sno 

IER=0 

CALL LORAN (ALAT,ALON,L,DEL,TD,ISYS,BETA,BLEN, BLEM, ISENT,ISW,IB, 


TIERS IT 31 el 


WRITE(9,20) I,H,M,S,ALAT(5),ALON(5),IT,IER 
FORMAT(13,1%,12,1X,12,1X,F5. 2.2%, Fl002, 23s ee 
F=ALAT(5)/3600. ODO 

LN=-ALON(5)/3600. ODO 
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WRITE(10,30) I,H,M,S,F,LN 
POnwet(13,1%,12,1X,12,1X,F5.2,2X,F10.7,2X,F12. 7) 
READ(8,40) F1,LNI 

FORMAT( 9X,F10.7,1X,F12. 7) 

DF=F-Fl 

DL=LN-LN1 

WRITE(6,*) DF,DL,LN1 

CONTINUE 

STOP 

END 


SUBROUTINE LORAN ( ALAT,ALON,L,DEL,TD,ISYS,BETA,BLEN, BLEM, 


*ISENT, ISW,1IB, IER, IT,IVEL) 


DATA DECLARATION 


IMPLICIT REAL*8(A-H,O-Z) 
INTEGER SOUTH, EAST ,WEST 


STORAGE ALLOCATION 


PEN onOn URI 2), 10( 2), BETACZ) ,BLEN(2) ,BLEM(2) , ISYS(2),HC2) ,R(2), 


*TVEL(2) , VELL(2) ,L(2) 


DIMENSION VEL(3) 

DIMENSION TEMP1(4) ,TEMP2(4) 

DIMENSION ALAT(5) ,ASON(5) 

DIMENSION A(9),BC9) 

Weems; 65/6137.0D0 , 6377397. 1550D0 , 6378293. 6450D0 


© 6378206. 4000D0 , 6378249. 1450D0, 6377276. 3450D0, 6378166. OOO00DO, 
*6378388. 0O000D0, 6378135. 0000D0/ 


Wat ee G256752,51420D0 ,6356078. 96Z28D0 ,6356617. 9376D0, 


Bae oso oU00N0.6356514. 8695D0,6356075. 4131D0,6356784. 2836D0, 
‘eegoo ble 946))0,6356750.5200D0/ 


DATA A1/24. 0305D0/ ,A2/-0. 40758D0/ ,A3/0. 00346776D0/ 
DATA B1/0. 510483D0/,B2/-0. 011402D0/,B3/0. 001760D0/ 
DATA ARC1/0. 484813681110D-5/ 

DATA LORANC/1/ 

DATA FAZ/0. ODO/,IK/0/,TEMP1/4*0. ODO/ , TEMP2/4*0. ODO/ 
DATA ZERO/0. ODO/ 

DATA NORTH,SOUTH,EAST,WEST/1HN,1HS, 1HE, 1HW/ 

DATA MINUS/1H-/ 

DATA VEL/299. 6929D0,299. 6911624D0, 299. 708D0/ 

DATA BAZ/0. ODO/,JER/0/ 

DATA ALATS/O. 0D0/ ,ALONS/0. ODO/,N/0/ 


INITIALIZE ERROR CODE 
IER=0 

TEST SWITCH FOR BASELINE COMPUTATION 
IF (IB) 5,150,5 


BASELINE COMPUTATION 


a Te a I | 


C2202 @ 


Be ge hs 


ag Fs ap 1 o> We 


10 
15 


20 


Ze 


30 


35 


40 


45 


50 
oe) 


60 


65 
70 


US 


80 


TEST FOR VALID VFLOCITY OF PROPAGATION CODE FOR BOTH CHAINS 


Oo) 0 Te) 2 

iF (IVEL(1). Go? GOunoms 
CONTINUE 

GO TO 20 

TER=1 

Conmonsz0 


DETERMINE VELOCITY OF PROPAGATION FOR BOTH CHAINS 


DO 25 I=1,2 
IVEL(I)=IVEL(1)+1 
KK=IVEL( 1) 
VELL( I )=VEL(KK) 
A4SA=VELL( 1) 
A4B=VELL( 2) 


TEST FOR CHANGE IN SPHEROID CODE, FOR VALID SPHEROID VODE 
AND DETERMINE APPROPRIATE SPHEROID COSTANTS 


DO 30 J=1,2 

IF (1(J). GE. 1 ANDO RCO). LE oS pmGOnnorso 
GO TO 55 

CONTINUE 

IF (L(1).EQ.L(2)) GO TO 35 
GOTO 55 

K=L(1) 

IF (IK-K) 40,60,40 

DO 50 I=1,9 

IF (I-K) 50,45,50 

A5=A(1) 

A6=B(1) 

A55=A5**A5 
A10=(A55-A6*A6)/A55 

IK=K 

GO TO 60 

CONTINUE 

IER=2 

GO TO 320 


TEST FOR CHANGE IN FIXED STATIONS 


DO 80 I=1,4 

IF (TEMP1(1)-ALAT(I)) 70,65,70 
IF (TEMP2(1)-ALON(1)) 70,80,70 
DO 75 J=1,4 

TEMP 1( J)=ALAT(J) 
TEMP2(J)=ALON( J) 

eo) aie) 35 

CONTINUE 

GO TO 145 


TEST FOR TRIAD OR TETRAD CONFIGURATION AND ACTIVATE 
APPROPRIATE SENTINEL 


56 
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mqaan 
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c7- C2 C2 


co €2C2 


85 


90 
2/5 


100 


TOS 


110 
lS 
120 


i> 


130 


35 


140 


145 


1,500 


A12=ALAT( 1) -ALAT(3)+ALON( 1) -ALON(3) 
A12=DABS(A12) 

A13=ALAT( 2) -ALAT( 4)+ALON( 2) -ALON(4) 
A13=DABS(A13) 

IF (A12-0.001D0) 90,90,105 

IF (A13-0.001D0) 95,95,100 

IER=3 

GO TO 320 

Al1=-1. ODO 

co nO 110 

A11=1. ODO 


TEST FOR VALID OPERATION CODE 


Peielowe) 120,125,115 
Hem@low=2) 125,125,120 
IER=4 

GO TO 320 


COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
MASTER 1 STATION TO SLAVE 1 STATION 


CALL SODIN (ALAT(1),ALON(1),ALAT(2) ,ALON(2),L,A35,FAZ, BAZ, JER) 
AA10=A35/AGA 

M=1 

I=] 

GO TO 295 

AA11=A47 


COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
MASTER 2 STATION TO SLAVE 2 STATION 


CALL SODIN (ALAT(3),ALON(3) ,ALAT(4) ,ALON(4) ,L,A35,FAZ, BAZ, JER) 
CC10=A35/A4B 

M=2 

I=2 

GO2fO 295 

CC11=A47 


COMPUTATION OF BASELINE DATA 
BETA(1)=AA10+AA11 
BETA(2)=CC10+CC11 
DG 140 I=1,2 
BLEN(I)=DEL(I) 
BLEM( I )=2. ODO*BETA(1I)+DEL(I) 
TURN SWITCH FOR BASELINE COMPUTATION OFF 


IB=0 
RETURN 


TEST SWITCH FOR FORWARD OR INVERSE COMPUTATION 
COmmemieles, 190), 1SW 


FORWARD COMPUTATION 


CC) GaGa Cc) 


Cla Ga c) ClCaiCa Ga 


C2 C7 2 02 


ee Ge le 
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160 


165 


170 


175 


180 


185 


COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
SHIP TO MASTER 1 STATION 


CALL SODIN (ALAT(5) ,ALON(5) , ALAT( 1) ,ALON(1),L,A35 ,FAZ, BAZ, JER) 
A76=A35/A4A 

A72=A76 

M=1 

[=7 

GOO 275 

A77=A47 

A73=A77 


COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
SHIP TO SLAVE 1 STATION 


CALL SODIN (ALAT(5),ALON(5),ALAT(2),ALON(2),L,A35,FAZ,BAZ, JER) 
A74=A35/A4A 

M=1 

I=8 

GO TO 295 

A75=A47 


COMPUTE GEODETIC DISTANCE AND SALT-WATER RETARDATION FROM 
SHIP TO SLAVE 2 STATION 


CALL SODIN (ALAT(5),ALON(5),ALAT(4) ,ALON(4) ,L,A35,FAZ,BAZ,JER) 
A70=A35/A4B 

M=2 

I=9 

GOO 295 

A71=A47 


COMPUTE GEODETIC AND SALT-WATER RETARDATION FROM 
SHIP TO MASTER 2 STATION 


IF (All) 185,175,175 

CALL SODIN (ALAT(5),ALON(5),ALAT(3),ALON(3),L,A35,FAZ, BAZ, JER) 
A72=A35 /A4B 

M=2 

ii 

GO TO 295 

A73=A47 


COMPUTE LORAN TIME DIFFERENCES OF FIX 
TDC 1)=AA10+AA11+A744+A75-A76-A77+DEL( 1) 
TD( 2)=CC10+CC11+A70+A71-A72-A734+DEL(2) 
RETURN 


INVERSE COMPUTATION 


SOR Eee POS TONSOE = Si LP 


C2 C2 °C? 


C2 G2 Ca" Go ey is as | La Ii ep as a 


Ca Ga eae.) 


a Us a I ee el ee 


190 


r9> 


200 


205 


210 


ZA 5 
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ALATDR=ALAT(5) 
ALONDR=ALON(5) 


TEST FOR VALID MICROSECOND READINGS 


DO 195 I=1,2 

IF (TD(1). LE. BLEM(I). AND. TD(1). GE. BLEN(I)) GO TO 195 
IER=5 

GO TO 320 

CONTINUE 


CONVERT SHIPS POSITION TO RADIANS 


TF (GESENT)>205;, 205,200 
ALATR=ALAT(5 )*ARC1 
ALONR=ALON(5)**ARC1 


INITIALIZE IETRATION COUNTER 


IT=0 
GO 10 215 


TEST ITERATION COUNTER 


ici P eit .30)°G0 TO 215 
TER=6 
GO 10° 320 


COMPUTE SINES AND COSINES OF SHIPS POSITION 


SINS=DSINCALATR) 
COSS=DCOS( ALATR) 


COMPUTE GEODETIC DISTANCE, SIN AND COSINE OF AZIMUTH, AND 
SALT-WATER RETARDATION FROM SHIP TO SLAVE 2 STATION 


CALL SODIN (ALAT(5),ALON(5),ALAT(4) ,ALON(4) ,L,A35,FAZ,BAZ, JER) 
M=2 

I=3 

GOTO 2580 

C1=A35 

D1=C1/A4B 

C2=DSIN(CFAZ) 

C3=DCOS( FAZ) 

C101=A47 


COMPUTE GEODETIC DISTANCE, SIN AND COSINE OF AZIMUTH, AND 
SALT-WATER RETARDATION FROM SHIP TO MASTER 2 STATION 


CALL SODIN (ALAT(5),ALON(5),ALAT(3),ALON(3),L,A35,FAZ,BAZ,JER) 
M=2 

I=4 

GO TO 280 

C4=A35 

D4=C4/A4B 

C5=DSIN(FAZ) 


C2) (aca ©) 


le a Se fe a> Je 
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250 


PLES 


240 


C6=DCOS(FAZ) 
C104=A47 


COMPUTE GEODETIC DISTANCE, SIN AND COSINE OF AZIMUTH, AND 
SALT-WATER RETARDATION FROM SHIP TO SLAVE 1 STATION 


CALL SODIN (ALAT(5) ,ALON(5) , ALAT(2) , ALON(2),L,435,FAZ, BAZ, JER) 
M=1 

I=5 

GO TO 280 

C7=A35 

E7=C7/A4A 

C8=DSIN(FAZ) 

C9=DCOS(FAZ) 

C107=A47 

IF (All) 245,235,235 


COMPUTE GEODETIC DISTANCE, SIN AND COSINE OF AZIMUTH, AND 
SALT-WATER RETARDATION FROM SHIP TO MASTER 1 STATION 
(FOR TETRAD CONFIGURATIONS ONLY) 


C10=C7 
C11=C8 
C12=C9 

CALL SODIN (ALAT(5),ALON(5) ,ALAT( 1) ,ALON( 1) ,L,A35 ,FAZ, BAZ, JER) 
M=1 

I=6 

GO TO 280 
C7=A35 
D7=C7/A4A 
C8=DSIN(FAZ) 
C9=DCOS( FAZ) 
C108=A47 


DETERMINATION OF TIME DIFFERENCE COSTANTS AND DIFFERENTIAL 
CORRECTIONS IN Y AND X 
(FOR TETRAD CINFIGURATIONS ONLY) 


C13=TD(2)-CC10-CC11-C101+C104-DEL(2) 
C17=C13*A4B 

C18=TD( 1)-AA10-AA11-C1074+C108-DEL(1) 
C22=C18**A4A 

C23=C1-C17 

C24=C4 

C25=C7+C22 

C26=C10 
C27=((C2-C5)*(C25-C26)*(C23-C24)*(C11-C8)) 
C29=((C2-C5)*(C12-C9 )*(C3-C6)*(C8-C11)) 
C30=C27/C29 
C28=(C23-C244C30*%(C3-C6))/(C5-C2) 

GO 10 250 


DETERMINATION OF TIME DIFFERENCE COSTANTS AND DIFFERENTIAL 


CORRECTIONS IN Y AND X 
(FOR TRIAD CINFIGURATIONS ONLY) 


60 
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275 
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C13=TD( 2) -CC10-CC11-C101+C104-DEL( 2) 
C17=C13*A4B 
C18=TD(1)-AA10-AA11-C107+C104-DEL( 1) 
C22=C18*A4A 

C23=C1-C17 

C24=C4 

C25=C7-C22 

C27=C2*(C25-C24 )4+C23*(C5-C8)+C8*C24-C5*C25 
C29=C2*(C6-C9)+C3*(C8-C5 )4+C5*C9 -C8*C6 
C30=C27/C29 
C28=(C23-C24+C30*(C3-C6))/(C5-C2) 


DETRMINATION OF DIFFERENTIAL CORRECTIONS IN LATITUDE AND 
LONGITUDE AND PRELIMINARY GEOGRAPHIC POSITION OF SHIP 


C100=1. ODO-A10*SINS*SINS 
©102—DSORT(C100) 

Se2—15)7/ 6102 
C31=C32*(( 1. ODO-A10)/C100) 
Gee —-650/001 

C34=( -C28/(C32*COSS)) 
ALATR=ALATR+C33 
ALONR=ALONR+C34 

Aieatees) —ALATR/ARC1 

ALON(5 )=ALONR/ARC1 
WRITE(13,*) ALAT(5),ALON(5),'ALAT ALON FM LORAN' 


TEST OF CONVERGENCE OF SOLUTION 


IF (All) 260,255,255 
H(1)=AA10+AA11+E7+C107-D7-C108+DEL(1) 
GO TO 265 
H(1)=AA10+AA11+E7+C107-D4-C104+DEL(1) 
H(2)=CC10+CC11+D1+C101-D4-C104+DEL( 2) 
fo) 270) 

RC I)=TD(1)-H(1) 
ANUM=DABS(R( 1) )+DABS(R(2)) 

IF (ANUM. LT. 0.0001D0) GO TO 275 
rn | 

GOMTO 210 

RETURN 


CONVERSION OF AZIMUTHS FROM NORTH TO SOUTH ORIENTATION 


ime AZ=643000. O00) 285,285,290 
FAZ=( FAZ+648000. 0D0)*ARC1 

Soto ee 5 

FAZ=( FAZ-648000. 0D0)*ARC1 


FORWARD AND INVERSE COMPUTATION 
OF SALT-WATER RETARDATION (FOR LORAN C ONLY) 
IF (ISYS(M). EQ. LORANC) GO TO 300 


A47=0. ODO 
Gort St5 


61 


moo moo 


Sp a a il 


300 A46=A35/1609. 349D0 

IF (A46-100.0D0) 305,310,310 
305 A47=(B1/A46)+B2+B3*A46 

GO TO 315 
310 A47=(A1/A46 )+A2+A3*A46 


BRANCHING FOR QUANTITIES TO BE USED IN SOLUTION 
315 GO TO (130,135,220, 225,230,240,160,165,170,180),I 
PRINT ROUTINE FOR INVALID DATA 


320 IF (IER. LE.4) GO TO 340 
CALL ANGLE (N,ALONDR, ISGN, LOND, LONM, ALONS) 
CALL ANGLE (N,ALATDR, ISGM, LATD, LATM,ALATS ) 
IF (ISGM. EQ. MINUS) GO TO 325 
KLA=NORTH 
GO TO 330 

325 KLA=SOUTH 

330 IF (ISGN. EQ. MINUS) GO TO 335 
KLO=WEST 
GO TO 340 

335 KLO=EAST 

340 GO TO (345,350,355, 360, 365,370), IER 

345 PRINT 375 
RETURN 

350 PRINT 380 
RETURN 

355 PRINT 385 
RETURN 

360 PRINT 390 
RETURN 

365 PRINT 395,LATD,LATM,KLA, LOND, LONM,KLO,TD(1),TD(2) 
RETURN 

370 PRINT 400,LATD,LATM,KLA, LOND, LONM,KLO,TD(1),TD(2) 
RETURN 


FORMAT STATEMENTS 


375 FORMAT(1H1,1X,' INVALID VELOCITY OF PROPAGATION CODE, DATA NOT PROC 
*ESSED' ) 

380 FORMAT(1H1,1X,' INVALID SPHEROID CODE,DATA NOT PROCESSED’ ) 

385 FORMAT(1H1,1X,' INVALID CONFIGURATION, DATA NOT PROCESSED’ ) 

390 FORMAT(1H1,1X,'INVALID OPERATION CODE, DATA NOT PROCESSED’ ) 

395 FORMAT(18X,213,A1,6X,14,13,A1,16X,F9. 2,13X%,F9. 2,35, INVALID RESDO 
*GS') 

400 FORMAT(18X,213,A1,6X,14,13,A1,16X,F9.2513x,F9.2,3X, INVALIDepReres 
TONE) 

END 


SUBROUTINE SODIN (ALAT1,ALON1,ALAT2,ALON2,L,S,FAZ,BAZ, IER) 
IMPLICIT REAL*8(A-H,0-Z) 


DIMENSTHONGS G2 E22 ht) 2) 
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DATA PI/3. 14159265359D0/ ,ARC1/0. 484813681110D-5/ 
DATA A,B/6378137. ODO, 6356752. 3142D0/ 


INITIALIZE ERROR CODE AND DETERMINE COSTANTS 


A14=1. OD0-(B/A) 
A56=A14*A14 
A50=1. OD0+A14+A56 
A51=A50-1. ODO 
A52=A56/2. ODO 
A53=-A51/2. ODO 
A54=A56/16. ODO 
A55=A56/8. ODO 
A57=A56%1. 250D0 
A58=A56/4. ODO 


TEST FOR IDENTICAL STATIONS 


IF(ALAT1-ALAT2) 45,35,45 
IF(ALON1-ALON2) 45,40,45 
IER=2 

GOerO 255 


COMPUTATION OF PARAMETRIC LATITUDE OF STATIONS 


X1=ALAT1*ARC1 
Y1=ALON1**ARC1 
X2=ALAT2**ARC1 
Y2=ALON2*ARC1 
X(1)=X1 
X(2)=X2 
BOrSs0>=1,2 
BETA=DATAN( (1. ODO-A14)*DSINC XC I) )/DCOSC X(T) )) 
EC 1I)=DSIN( BETA) 
F(I)=DCOS( BETA) 


COMPUTATION OF SPHERICAL DIFFERENCE OF LONGITUDES 


A59=-Y2 

A60=-Y1 

C35=A59-A60 
C36=DABS(C35) 
meG2o-F 1) )60555,55 
A16=2. ODO*PI-C36 
COeTO 65 

A16=C36 

ie@ele) 75,/0,75 
Al6=0. 5D-7 


COMPUTATION OF GEODETIC DISTANCE 


A17=DSIN(A16) 
A18=DCOS(A16) 
A19=E(1)**E(2) 
A20=F(1)*F(2) 
A21=A19+A20%*A18 
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A80=A17*F(2) 

A81=A80*A80 

A82=E(2)*F(1)-E(1)*F(2)*A18 

A83=A82**A82 

A84=A81+A83 

A22=DSQRT(A84) 

A23=(A20*A17)/A22 

A24=1. OD0-A23*A23 

A90=A22/A21 

A25=DATAN(A90) 

IF(A25) 80,85,85 

A25=(PI/2. 0D0)+A25 

A26=A25*A25 

A27=1. 0D0/A22 

A91=A21/A22 

A29=A24*A24 

A30=( A50%*A25 )+A19%*(A51*A22-A52*A26*A27) 
A31=A247*( A53*A25+A53%A22°A214+A52*A26*A9 1) 
A32=A19*A19%( -A52*A21%A22) 

A33=A29*( A54%A25+A54°A22*A2 1 -A52*A26%A9 1-A55*A22*A2 13) 
A34=A19*A24%( A52*A26%A27+A52*A22*A21%*A21) 
S=(A30+A31+A32+A33+A34)*B 


COMPUTATION OF DIFFERENCE OF LONGITUDE ON THE REDUCED SPHERE 


A36=(A51*A25 )+A19%*( -A52*A22-A56*A26*A27) 
A37=A24*( -A57*A254+A58*A22*A21+A56*A26*A9 1 ) 
A38=(A36+A37 )*A23+A16 


COMPUTATION OF GEODETIC FORWARD AZIMUTH 


A39=DSIN(A38) 
A4O=DCOS(A38) 
AG1=(E(2)*F(1)-A40*E(1)*F(2))/(A39*F(2)) 
IF(A41) 95,90,95 
A41=0. 5D0-7 

A42=1. 0D0/A41 
A43=DATAN( A42) 

IF(C35) 120,100,100 
TF(GS5"P9)) 105,115 105 
IF(A41) 110,140,140 
A43=PI+A43 

GO TO 140 

IF(A41) 130,135,135 
IF(C35+PI) 105,105,125 
IF(A41) 130,135,135 
A43=PI-A43 

GO TO 140 

A43=2. ODO*PI-A43 
A43=A43+PI1 

A43=A43-2. ODO*PI 
IF(A43) 145,150,150 
A43=A43+2. ODO*PI 
FAZ=A43/ARC1 
FAZ=FAZ+648000. ODO 
IF(FAZ-1296000. ODO) 160,155,155 
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FAZ=FAZ-1296000. ODO 
COMPUTATION OF BACK AZIMUTH NOT INCLUDED IN THIS VERSION 


BAZ=0. ODO 
RETURN 


PRINT ROUTINE FOR INVALID DATA 


Pela 270 

CALL ANGLE(3,ALAT1,ISGN1,LAT1D, LAT1M, ALATIS) 
CALL ANGLE(3,ALON1, ISGN2,LON1D, LON1M,ALON1S ) 
CALL ANGLE(3,ALAT2,ISGN3, LAT2D, LAT2M, ALAT2S ) 
CALL ANGLE(3,ALON2,ISGN4, LON2D, LON2M, ALON2S ) 
PRINT 280,ISGN1,LAT1D, LAT1M,ALAT1S,ISGN2,LON1D,LON1M,ALON1S,ISGN3, 


*LAT2D, LAT2M, ALAT2S , ISGN4 , LON2D , LON2M, ALON2S 


RETURN 


FORMAT STATEMENTS 

FORMAT( 1H1,51X,'SUBROUTINE SODIN INVALID DATA'/) 

FORMAT( 1X,4( 2X, Al, 213,F7. 3),2X,12,2X,' IDENTICAL STATIONS' ) 
END 


SUBROUTINE ANGLE (N,ARGS,ISIGN, IDEG,MIN,SEC) 


IMPLICIT REAL*8(A-H,0-Z) 
DIMENSION ASEC(6) 
Dat AP ASEG/ 594500 ,59. 945D0,59. 9945D0,59. 99945D0, 59. 999945D0, 


*59.9999945D0/ 


DATA MINUS/1H-/,IBLANK/1H / 
ANGLE CONVERSION 


IF (ARGS) 10,5,5 
ISIGN=IBLANK 
BRGS=ARGS 

GO TO 15 
ISIGN=MINUS 
BRGS=-ARGS 
IDEG=BRGS/3600. ODO 
ARGT=IDEG* 3600 
AMIN=BRGS-ARGT 
MIN=AMIN/60. ODO 
AMINT=MIN*60 
SEC=AMIN-AMINT 
M=N+1 

IF (SEC-ASEC(M)) 30,20,20 
SEC=0. ODO 
MIN=MIN+1 

IF (MIN-60) 30,25,25 
MIN=O 

IDEG=IDEG+1 

RETURN 

END 


APPENDIX D. TABLE FOR STANDARD DEVIATIONS AND 
COVARIANCES OF VELOCITY COMPONENTS FOR MINI RANGER 


TIME IN- 
Oy Ou Ov,u 
TERVAL 
3 MIN 1.3 cm/sec 0.1234 cin} sec? 


6 MIN 0.64 cm/sec 0.64 cm/sec 0.0336 Ge see- 


0.32 cimsec 32 einsce 0.0116 crm?/ sec? 
0.19 cmysec Gri orem see 10-' x 0.04 8G.em- sec: 
10-3 x 0.98cr/ sec! 10-* x 0.97 cm] sec OO Gr sees 
10-3 x 0.67cn/ sec | 10-3 x 0.66cni] seen 00 0s 0c77 eec- 
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APPENDIX E. AVERAGE CURRENT PROFILES 
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Figure 17. AVERAGE CURRENT PROFILE FROM POINT 1 TO POINT 3 
USING NAVIGATION DATA FROM MINI RANGER 
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Figure 18. AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 2 


USING NAVIGATION DATA FROM MINI RANGER 
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AVERAGE CURRENT PROFILE FROM POINT 2 TO 
USING NAVIGATION DATA FROM MINI RANGER 
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AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 4 
USING NAVIGATION DATA FROM MINI RANGER 
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VELOCITY (cn/s) 
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Figure 21. 
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AVERAGE CURRENT PROFILE FROM POINT 1 TO POINT 3 
USING NAVIGATION DATA FROM LORAN (TD) 
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AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 2 


USING NAVIGATION DATA FROM LORAN (TD) 
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VELOCITY (cn/s) 
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AVERAGE CURRENT PROFILE FROM POINT 2 TO POINT 3 


USING NAVIGATION DATA FROM LORAN (TD) 
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Figure 24. 
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AVERAGE CURRENT PROFILE FROM POINT 3 TO 
USING NAVIGATION DATA FROM LORAN (TD) 
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Figure 25. AVERAGE CURRENT PROFILE FROM POINT 1 TO POINT 3 
USING NAVIGATION DATA FROM LORAN (DISPLAY) 
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Figure 26. AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 2 
USING NAVIGATION DATA FROM LORAN (DISPLAY) 
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AVERAGE CURRENT PROFILE FROM POINT 2 TO POINT 3 
USING NAVIGATION DATA FROM LORAN (DISPLAY) 
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Figure 28. AVERAGE CURRENT PROFILE FROM POINT 3 TO POINT 4 


USING NAVIGATION DATA FROM LORAN (DISPLAY) 


78 


APPENDIX F. CURRENT PROFILES FOR SPECIFIC TIMES 
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Figure 29. CURRENT PROFILE FROM POINT | TO POINT 2 USING 
NAVIGATION DATA FROM MINI RANGER FOR A SPECIFIC 
TIME (14 33 11) 
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Figure 30. CURRENT PROFILE FROM POINT 2 TO POINT 1 USING 
NAVIGATION DATA FROM MINI RANGER FOR A SPECIFIC 


TIME (15 06 10) 
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Figure 31. CURRENT PROFILE FROM POINT 1 TO POINT 2 USING 
NAVIGATION DATA FROM LORAN (TD) FOR A SPECIFIC 
TIME (14 33 11) 
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Figure 32. CURRENT PROFILE FROM POINT 2 TO POINT 1 USING 
NAVIGATION DATA FROM LORAN (TD) FOR A SPECIFIC 


TIME (15 06 10) 
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Figure 33. CURRENT PROFILE FROM POINT 1 TO POINT 2 USING 
NAVIGATION DATA FROM LORAN (DISPLAY) FOR A SPE- 
CIFIC TIME (14 33 11) 
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Figure 34. CURRENT PROFILE FROM POINT 2 TO POINT 1 USING 
NAVIGATION DATA FROM LORAN (DISPLAY) FOR A SPEI- 
CIFIC TIME (15 06 190) 
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Compare at sea posi- 
tion using Mini-Ranger, 
Loran C (Internav) in 
the context of measur- 
ing current velocity 
with a shipboard ADCP 
(Acoustic Dopller [i.e. 


Doppler] Current Pro- 
filer). 





